Many body effects in the excitation spectrum of a defect in SiC 
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We show that electron correlations control the photophysics of defects in SiC through both renor- 
malization of the quasiparticle bandstructure and exciton effects. We consider the carbon vacancy, 
which is a well-identified defect with two possible excitation channels that involve conduction and 
valence band states. Corrections to the Kohn-Sham ionization levels are found to strongly depend 
on the occupation of the defect state. Excitonic effects introduce a red shift of 0.23 eV. The analysis 
unambigiously re-assigns excitation mechanism at the thresholds in photo-induced paramagnetic 
resonance measurements [J. Dashdorj et al, J. Appl. Phys. 104, 113707 (2008)]. 
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Optical spectra of semiconductors not only contain the 
excitations of the perfect bulk crystal but also that of its 
omnipresent imperfections, such as point defects. Pho- 
toluminiscence plays a fundamental role in the experi- 
mental characterization and identification of defects 
Despite the success of density functional theory (DFT) 
in explaining the physics of defects in covalent semicon- 
ductors [2], this common approach has severe shortcom- 
ings: the description of the bonding of the bulk crys- 
tal may be insufficient within common approximations 
for exchange and correlation, hke the local spin density 
approximation (LSDA), the position of defect levels in 
the band gap is affected by the well-known Kohn-Sham 
band gap error, and excitations cannot be assessed rig- 
orously. Many-body perturbation theory allows to re- 
solve these issues [3] , however, applications to defects are 
scarce ^ due to computational complexity. The quan- 
titative prediction of absorption spectra of a wide-class 
of systems from insulators to surfaces, to nanotubes and 
polymers became tractable only recently by solving the 
Bethe-Salpeter equations (BSE) and accounting for the 
Go Wo self-energy OE]. 

In electron paramagnetic resonance (EPR) experi- 
ments under illumination [BHEj, the positively charged 
carbon vacancy (Vj) in the 4H polytype of silicon car- 
bide (4H-SiC) was extensively investigated and excitation 
thresholds were assigned to the defect ionization levels 
with respect to the neutral defect. The identification of 
this defect as an EPR-center was established on the basis 
of the DFT-LSDA ground state calculations jlllU]. Yet, 
for the accurate interpretation of the photo-EPR data 
one has to consider excitation transitions taking into ac- 
count at least two competing excitation channels: tran- 
sitions to the neutral state (V,'.) as well as to the dou- 
bly positive defect (V^^). To assess the values of the 
ionization levels one has to evaluate the contribution of 



excitonic effects that is unknown so far. Knowledge of 
the ionization levels is pivotal for an understanding of 
the carrier compensation in high purity semi-insulating 
SiC [TT] , in which vacancies are abundant compensation 
centers, and for the dopant diffusion and defect kinetics 

m- 

Defect excitations combine the many-body effects in 
both extended and localized states. The insertion of elec- 
trons into defect states invokes on-site correlation due to 
electron-electron repulsion expressed by the Hubbard U. 
Defect-to-band transitions involve two particles, the elec- 
tron and the hole (one in the localized state and the other 
in the extended state). To self-energy effects this adds 
excitonic correlation energy of yet unknown size. 

The aim of this letter is to investigate many-body ef- 
fects in the absorption spectrum of a well-identified de- 
fect: a carbon vacancy in 4H-SiC. We show that those 
corrections strongly depend on the occupation of the de- 
fect levels, indicating a Hubbard-?/ which is substan- 
tially larger than predicted by DFT-LSDA (1.09 vs. 
0.79 eV). Furthermore excitonic effects introduce a size- 
able red shift (~ 0.23 eV) varying with the defect charge 
state. More importantly, the calculated transition en- 
ergy thresholds for the two competing channels indicate 
that the current interpretation of photo-EPR spectra [B]- 
IH] based on the assumption of a photoionization via the 
neutral defect has to be corrected. According to our 
calculations the ionization to V^"*" is the only excitation 
mechanism at photon energies below 2.3 eV that is re- 
sponsible for the experimental findings. 

The defect absorption spectra and quasiparticle (QP) 
energies were investigated with the YAMBO-code |13j . 
The defect in 4H-SiC was represented by a supercell with 
288 lattice sites [14]. Among the two inequivalent va- 
cancy sites in 4H-SiC [9] we focus on the experimentally 
investigated cubic site. In the GqWo and BSE calcu- 
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FIG. 1: Carbon vacancy Vj in 4H-SiC (color on-line): (a) ge- 
ometry - the spin density is indicated, (b) defect levels within 
the band gap. In the ground state, the a-level is occupied with 
the paramagnetic electron. Arrows indicate the transitions 
from the valence band (V) (in orange) and to the conduction 
band (C) (in green), (c) sketch of the total energy i?tot(Q) 
of Vj and Vc vs. the atomic coordinates; the QP-energies 
(-|-|0)v and (0|-|-)v and the Franck-Condon shift Eyc of V% 
are indicated (cf. text). 



TABLE I: QP-energies {q\q + l)v and {q\q - l)v, Franck- 
Condon shifts Eyc and thermal ionization levels {q\q — 1) for 
the relevant charge states of Vc in eV. Values are listed as 
obtained within Go Wo and DFT. First two columns: initial 
charge state with the occupation of its defect levels. Efc is 
taken from DFT calculation after electron insertion (cf. text). 



(gk + i)v (gk-i)v Eve (gk-i) 

GoWo DFT GqWo DFT GoWo DFT 

1.83 1.85 0.29 1.54 1.56 

1.27 1.24 2.36 2.03 0.83 1.53 1.20 

0.72 2.80 2.28 0.21 2.59 2.07 

1.39 2.65 2.16 0.41 2.24 1.76 



lations spin was taken into account |15| . First we turn 
to the fundamental band gap of 3C and 4H-SiC. The 
QP-corrections using the DFT-lattice constant amount 
to 0.91 eV and 0.93 eV respectively, in agreement with 
earlier GoWo-calculations 16J. Thereby, for the gap of 
4H-SiC, we obtain 3.11 eV (Exp. 3.265eV [Hj). 

Optical absorbtion, as indicated in Fig. [T] involve 
the defect levels in the band gap, a non-degenerate 
(a-level) and a two-fold degenerated level (e- level). A 
(pseudo) Jahn- Teller effect [S] lifts the degeneracy for 
all charge states except for V^"*", however, for Vj the 
effect is neglegible. Electronic transitions occur predom- 
inantly at the fixed geometry of the initial state (Franck- 
Condon principle). The excitation threshold for tran- 
sitions between the valence band and the defect states 
is determined by the QP-energies of the lowest unoccu- 
pied defect state (LUMO) Elumo ^'^'^ valence band 
edge Ey and by the electron- hole interaction i^ohj i-e. 
^LUMO ~ ~ Ech- For transitions from the highest oc- 
cupied defect state (HOMO) to the conduction band with 



Ionization - We now discuss the QP-contribution to 
the threshold in comparision with the DFT values to dis- 
cern the many-body corrections within the GqWo- The 
addition of an electron to the LUMO of the defect in the 
charge state q (V^) yields V^~^. Using E^ as an en- 
ergy reference for the electron, the energy required for 
the addition is described by the QP-energy difference 

-LUMO 



i?Vj which we denote {q\q— l)v Correspond 



ingly, the removal of an electron from the HOMO of Vq, 
which yields Vq^^, is associated with the energy gain 
^HOMO ~ denoted by {q\q + l)y. In DFT, the val- 
ues are obtained from the total energies. The results are 
listed in Table [ij Consider first V^"^ and Vj: for the 
addition of an electron to the unoccupied a-level of V^"^ , 
(2-)-|-|-)v obtained within GqWo and DFT agree to within 
20meV. Also, for the electron removal from the a-level 
of Vj, (+|2-h)v obtained with GqWo matches the DFT- 
result. Hence, the addition (removal) of an unpaired 
electron to (from) the a-level is already well described 
within DFT. Here the QP-correction to the Kohn-Sham 
level of 0.2 eV (insertion) and -0.23 eV (removal) reflects 
the screening by valence band electrons and defect res- 
onances. However, for adding a second electron to the 
a-level exchange and correlation effects beyond the de- 
scription of the DFT considerably change the mutual re- 
pulsion of the two electrons: The GqWo value for (+|0)v 
is by 0.30 eV larger than the DFT result. The effective 
Hubbard interaction U = (+|0)v - (+|2+)v (cf. [18] ) 
among the two electrons in the QP-treatment amounts 
to [/ = 1.09 eV instead of the DFT value U = 0.79 eV. 
For the remaining ionization levels, similar changes in 
the electron-electron repulsion due to the inclusion of ex- 
act exchange and the screended Coulomb-interaction in 
GqWo beyond the DFT description are found. 

The thermal ionization level {q\q— 1) determines in 
thermal equilibrium the positions of the Fermi level at 
which the charge state changes from qio q — 1. The level 
{q\q~ 1) differs from {q\q — l)y by the Franck-Condon 
shift, which is the energy gain due to the relaxation 
of the vacancy after the electron insertion (cf. Fig. [ij;). 
The relaxation leads to an effective electron- electron re- 
pulsion Ues smaller than U. To assess {q\q— 1) within 
the GgWo-approach, we approximate the Franck-Condon 
shift E-pc by the DFT values. The results are listed in 
Table |l] For Vj , DFT predicts a negative value for UcS 
p] . Within DFT we obtain C/es = -0.36 eV. This would 
imply that the charge state is unstable and in contra- 
diction to the findings [5J [7] should hardly be observable 
by EPR in thermal equilibrium. However, due to the QP- 
corrections, the levels (2-1- |-|-) and (+|0) become almost 
equal (i.e. UcS ~ 0), thus resolving the contradiction. 
Similarly, a negative-?/ effect is found for V^. Here, af- 
ter the QP-correction, Ucs remains negative with a value 
of -0.35 eV. 
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FIG. 2: (color online) Absorption spectra of Vc within 
GoWo-BSE: (a) for Vj in comparision with the GqWo-RPA 
spectrum and the joint densisty of states (jDOS) of the tran- 
sistions with the a-level, (b) for Vc and V^^ as calculated 
with spin-polarization and without (for the latter). The on- 
set of excitation channels is indicated. A broadening of 0.1 eV 
was used for all spectra. 



TABLE II: Ionization (Go Wo) and excitation (BSE) thresh- 
olds as well as experimental data. Also the excitation energy 
of transition between the a and e level are given. Values are 
in eV. 



V2+ 



Go Wo BSE £;eh 



a - 
V- 
V- 
a - 



■ e 
a 
e 

■ C 



1.83 
2.14 



1.65 0.18 
1.96 0.18 



GoWo BSE Boh 
1.09 0.89 0.20 
2.39 2.16 0.23 
2.36 2.13 0.23 
1.84 1.61 0.23 



Go Wo BSE E,i, 
1.92 1.68 0.24 

2.80 2.56 0.21 
2.23 2.02 0.21 



Experimental thresholds 






quenching of V J 


1.47" 


1.6' 


restoring of Vj 


1.81" 


1.9' 



"Ref. [3 
''Ref. [8] 



Optical properties ~ Photo-EPR experiements [BHl] 
were conducted with Vq as the initial or the final state. 
The electronic transitions, that quench the paramagnetic 
state and include the ionization of Vj into Vq"^ or Vq, 
are shown in Fig. [ij (i) excitation from the a-level to the 
conduction band (a — >■ C), or (ii) from the valence band 
to the defect levels a or e {V ^ a or V ~^ e). Restoring 
the EPR-signal of starts from V^"^ or as initial 
states by exciting an electron from the valence band to 
the a-level (V^"^) or from the a-level to the conduction 
band (Vq) and implies the subsequent dissociation of the 



electron-hole pair. We calculated absorption spectra for 
Y'q~ , Vq and, to unravel the relevance of competing 
channels and to assess the electron-hole coupling. For 
the latter purpose we compared the spectra with the in- 
clusion of the electron-hole interaction via the BSE and 
without it, in the random phase approximation (RPA). 
In the investigated energy range only absorbtion via de- 
fect states contributes to spectra due to the indirect band 
gap. The spectra are shown in Fig. [2][T1]. In Table |TT) 
we list the ionization and excitation thresholds. We also 
calculated RPA spectra, shown in Fig. |3j for the denser 
4x4x4 k-point mesh [20] that are converged with respect 
to the density of extended states. Since such calculations 
are not feasible for the BSE, corresponding BSE spec- 
tra were extrapolated by a convolution ansatz. Spectral 
broadening due to non-vertical transitions is accounted 
for by applying a Huang-Rhys factor for the coordinate 
that connects the geometry of the ground state and the 
ionized state (broading amounts to ^0.26 eV). 

First we address the spectra of obtained for the spe- 
cial k-point (Fig. [2^). Excitonic effects lead to an almost 
rigid red shift of the RPA-spectrum by E^h = 0.23 eV. 
The same effect is found also for Yq, and Vq albeit 
with different values of E^h (cf. Table |lT| . The structure 
of the RPA-spectrum is dominated by band structure ef- 
fects. This is demonstrated by comparison with the joint 
density of states (jDOS) of the transitions involving the 
a-level, which essentially reproduces the RPA-spectrum. 
The spectral broadening by Huang-Rys factors, however, 
reduces these features as seen in Fig. [3j Deviations be- 
tween the BSE and RPA spectra, besides the red-shift, 
stem from bright and dark exciton solutions of the BSE 
due to the spin-dependent electron-hole exchange cou- 
pling. The reduction of absorption due to dark exciton 
states is, however, not obtained in a spin-averaged treat- 
ment. This is seen in Fig. [2|d, where we compare the 
BSE-spectra of V^"^ as obtained in the two ways. 

The vertical excitation a — >■ e between the defect levels 
occur at 0.89 and 1.68 eV for and V'^ well below the 
onset of transitions involving extended states. An im- 
portant result is that the transitions to the conduction 
band set in at energies (vertical transitions at 1.61 and 
2.02 eV for and V° , respectively) below the threshold 
for excitation from the valence band (2.13 eV for and 
2.56 eV for V^). 

From recent time-resolved photo-EPR experiments [2] 
cross sections for quenching and restoring of the EPR sig- 
nal of are available with thresholds at 1.6 and 1.9 eV 
respectively. Steady state photo-EPR [7] yields values of 
1.47 and 1.81 eV, respectively (cf. Table |ll|. The calcu- 
lated onset for excitation a — > C at 1.61 eV agrees well 
with the experiments. A rise of the quenching cross sec- 
tion at 2.4 eV [8] should be related to our finding for the 
onset of the transitions V a at 2.16 eV and direct ion- 
ization above 2.4 eV that yields V^. For restoring V^, 
however, the excitation of an electron-hole pair with its 
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FIG. 3: (color online) Absorption spectra of Vj and Vq^. 
GoWo-RPA and extrapolated GqWo-BSE spectra for high It- 
point density. The effect of Franck-Condon shifts are shown 
(indicated by the extension HR). A broadening of 0.01 eV was 
used for the raw GoWo-RPA spectra. 

subsequent thermal dissociation or direct ionization is re- 
quired with thresholds of 1.65 and 1.83 eV, respectively. 
The comparison with experiment indicates that direct 
ionization is the relevant process. The present analysis 
of competing excitations suggests a re-interpretation of 
the photo-EPR spectra [BHH] in terms of a photoioniza- 
tion of via V^"^ instead of Yq as considered there. 
Consequently, the ionization levels deduced from the ex- 
periments have to be reassigned correspondingly and cor- 
rected for the exciton binding energy. 

In summary we conducted GW-I-BSE full fledge calcu- 
lations of the photophysics of the carbon vacancy in 4H- 
SiC. We demonstrate the importance of adding correla- 
tion effects at that level, in spite of the fact that the struc- 
ture is well described at the DFT-level. Charge state- 
dependent corrections include the ionization levels ( up 
to 0.5 eV), the electron-electron repulsion U (~ 0.3 eV) 
and the electron- hole attraction (~ 0.23 eV). This has 
implications for the analysis of defects in materials, such 
as oxides, where the correct despription of the bonding 
is an issue. We showed that ionization of the defect into 
Vq^ and Yq compete. However, for the range of pho- 
ton energies in which photo-EPR experiment identified 
thresholds only the former ionization channel is active. 
Excitation spectra are red shifted by strong excitonic ef- 
fects. This unambigously rectifies the earlier experimen- 
tal assignment of transitions and the assessment of ion- 
ization levels. These are pivotal for the defect diffusion 
and for unraveling the mechanisms carrier compensation 
in semi-insulating SiC and hence for applications. 
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